Sub-diffraction limited quantum imaging of a living cell 



Michael A. Taylor, 1 ' 2 Jiri Janousek, 3 Vincent Daria 4 , Joachim Knittel 1 , 
Boris Hage 3 , Hans-A. Bachor 3 ' 4 , and Warwick P. Bowen 2 * 

department of Physics, University of Queensland, St Lucia, Queensland 4072, Australia 
2 Centre for Engineered Quantum Systems, University of Queensland, St Lucia, 
Queensland 4072, Australia 
3 Department of Quantum Science, Australian National University, Canberra, 

ACT 0200, Australia 

4 John Curtin School of Medical Research, Australian National University, Canberra, 

ACT 0200, Australia 

*To whom correspondence should be addressed; E-mail: wbowen@physics.uq.edu.au 



Abstract 

Quantum techniques allow the resolution constraints of classical imaging to be overcome, and are expected 
to have important applications in biology. We report the first demonstration of sub-diffraction limited quantum 
imaging in biology. Naturally occurring lipid granules of approximately 300 nm diameter are used to image 
the local mechanical properties of the cellular cytoplasm, with spatial resolution enabled by thermal diffusion. 
Spatial structure is resolved at length scales down to 10 nm. Our results confirm the longstanding prediction 
that use of quantum correlated light can enhance spatial resolution in biology, allowing a 14% enhancement 
over that achievable with coherent light. Combined with state-of-the-art sources of quantum correlated light, 
this technique provides a path towards an order of magnitude improvement in resolution over similar classical 
imaging techniques. 

Quantum resources play an increasingly important role in precision measurements, and are applied in widely 
varying fields such as astrophysics [1 1, field sensing J2] [3] 0), and subcellular microscopy J5] |6] [8] . They 
both allow classical constraints such as the shot-noise limit or diffraction limit to be overcome, and enable new 
engineered nanoprobes for sensing. Such advances are important in biological imaging, where any improvement in 
imaging technology can reveal new levels of cellular complexity. Since sub-cellular structures can have nanometer 
size scales, spatial resolution surpassing the diffraction limit is particularly beneficial. The realization of sub- 
diffraction limited quantum imaging in biology is therefore a longstanding and important challenge. 

Here we demonstrate sub-diffraction limited biological imaging with quantum enhanced resolution. To achieve 
this result, a new approach to quantum imaging is developed which applies squeezed states of light within a 
classical sub-diffraction limited imaging technique. Together this allows cellular structure to be imaged at length 
scales down to 10 nm. For fixed optical power, squeezed light provides 14% enhancement over the resolution 
possible with classical imaging, demonstrating for the first time that non-classical light can improve resolution 
in biology. Alternatively, a 42% reduction in optical power is possible without compromising resolution, with 
consequent reduction in damage |j9j [TOl and photochemical disruption of cellular processes (TTOl ITTI . 

This first quantum image of the structure within a living cell complements other recent applications of quantum 
imaging in biology, with Nitrogen Vacancy (NV) nanodiamonds allowing thermal [7 1 and magnetic imaging |[6] 
fT2l n~3|| . and quantum correlated light allowing dispersion compensation 1 14]. Prior to the results presented here, 
all biological applications of quantum imaging have been diffraction limited, including NV based cellular imaging 
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which has relied on optically resolvable arrays of stationary NV probes ||6]|7][T2][T3], confining it to the study of 
relatively large cellular structures and organelles. The resolution achieved here is over an order of magnitude finer 
than in previous biological quantum imaging experiments, providing the possibility to observe important nanoscale 
cellular structures such as membranes, actin networks, and individual proteins. Since the approach is in principle 
transferrable to NV based imaging, it could also open the door to simultaneous sub-diffraction limited imaging of 
viscoelasticity, temperature and magnetic fields. 

In the imaging technique introduced here, the local structure of the cell is interrogated non-invasively via its 
effect on the thermal motion of naturally occurring nanoparticles as they diffuse around the cell. In a purely viscous 
medium such as water, these particles would follow well-known diffusive Brownian motion. By contrast, the cy- 
toplasm in a cell is structured and crowded with proteins and other macromolecules which confine particle motion 
and slow the thermal motion 1151 . The resulting subdiffusive motion characterizes the viscoelasticity of the cyto- 
plasm |[T6l[T7l[T9l[T8l . Thermal motion is crucial to the operation of the cell, mediating many important processes. 
For instance, subdiffusion slows the initial approach of chemical reactants to one another but also decreases the 
likelihood that they will pass without interacting, which improves the efficiency of the reaction [ 19, 20. 21]. Subd- 
iffusion also has a strong effect on enzyme reactions [22, 23], and the formation of spatiotemporal patterns ||24]|251 . 
The technology introduced here allows improved observation of spatial variations in subdiffusive motion, and is 
thereby an enabling step towards a better understanding of the nanoscale structure within a cell. Although the 
biological implications of such structure are not yet fully understood, it has been shown that the optimal diffusive 
regime is different for storage, transport, and chemical reactions lEUl . As such, nanoscale structural variations 
could potentially reflect the function of small regions of the cytoplasm. 

The experiment builds upon the recently developed quantum enhanced nanoparticle tracking technique de- 
scribed in Ref. [8|. In that work, temporal fluctuations in cellular viscoelasticity were monitored with quantum 
enhanced precision via optical centroid measurements [26. 27] . However, the possibility to use diffusion to achieve 
spatial resolution, and thereby image the cell, was not examined. As a result fluctuations caused by time-dynamic 
processes within the cell could not be distinguished from those due to nanoscale structure that the particle is ex- 
posed to as it diffuses around the cell. Nor was it possible to retrieve any quantitative information about the spatial 
structure. These limitations are resolved here by tracking nanoparticles with sub-nm precision as they explore 
an extended region of the cell, and analyzing spatial variations in the subdiffusive motion. The experimental 
setup (shown in Fig. [T]\) features several important modifications from a conventional optical particle tracking 
experiment, facilitating the use of amplitude squeezed light to enhance measurement sensitivity. Because most 
biological processes occur at Hz-kHz frequencies, where classical noise sources constrain the possibility of gen- 
erating squeezing l28l . a novel optical lock-in technique [29 1 is used to evade low frequency noise and allow 
quantum enhancement at these frequencies. This is combined with dark-field illumination ll30l[3T1 to remove un- 
wanted light from the measured optical field. Additionally, self-homodyne measurement is used, which allows 
the local oscillator itself to be squeezed, such that the squeezed field perfectly overlaps with the local oscillator 
at detection even without prior knowledge of the spatial mode shape after propagation through the living cell and 
high numerical aperture lenses. 

Saccharomyces cerevisiae yeast cells were immobilized with an optical trap, and the motion of lipid granules 
within the cellular cytoplasm was tracked with either squeezed or coherent light J8]. The high frequency thermal 
motion of a lipid granule then reveals the viscoelastic properties of the surrounding cytoplasm, while for sufficiently 
long measurements its slow thermal drift provides spatial resolution by bringing it into contact with different 
parts of the cellular cytoplasm. Thus, spatial inhomogeneity in the viscoelasticity can be quantified by a single 
continuous measurement of the lipid granule position. In our experiment, the particle position x(t) was measured 
by combining scattered light from the sample with a local oscillator field which was spatially shaped such that 
direct measurement of the total power yielded the particle position (Fig.[TJj). In the same manner as in Ref. J8j , 
the local viscoelasticity of the cytoplasm could be determined from the mean squared displacement (MSD) of the 
particle after a delay r, 



with an example shown in Fig.[Tfl. The corresponding spectra (Fig.QjT inset) shows that squeezed light lowers the 
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noise floor of the measurement by 2.4 dB. For short delays, the MSD is dominated by thermal motion and has the 
form 

(Ax 2 (r)) =2Dr a , (2) 

where the diffusive parameter a carries information about mechanical properties of the surrounding medium l32l 
[33l . a is determined for a set of data by fitting the MSD at short delays to Eq. [2] When a = 1, the motion 
is diffusive, which is indicative of a random walk type of motion, whereas confinement of the particle causes 
subdiffusive motion (0 < a < 1). In our experiments, 100 ms of data was sufficient to precisely determine a 
Subdiffusive motion is an indicator that the cellular cytoplasm exhibits both viscosity and elasticity lfT6l[T7l[T9l[T8l . 
since to constrict motion the cytoplasm must store mechanical energy. Consequently, the measured values of a 
allowed temporal variations in the cellular viscoelasticity to be characterized with 10 Hz bandwidth. 

As lipid particles undergo three dimensional (3D) thermal motion, they are exposed to different parts of the 
cell (Fig.QJ} inset). Since our experiments provide both a and the mean position of the particle along the x axis, 
it is possible to construct ID images of a which follow the projection of the 3D trajectory onto the x axis (see 
Fig. |2), and map out the physical structure within the cellular cytoplasm. To construct such images, a series of 
experiments were performed in which the nanoparticle motion was tracked as it diffused through an extended 
region of the cell. The data from each experiment was separated into 100 ms segments, with both a and the mean 
position along the x axis determined for each segment. As the particle diffused, an image of a could thereby 
be created as a function of x, as shown in Fig. |2] The directly obtained images exhibited substantial noise both 
from the measurement process and due to the unknown trajectory of the particle in the y and z directions. To 
identify statistically resolvable features, the running mean and standard error of a were calculated along the x axis, 
with spatial resolution given by the size of the averaging window. As shown in the supplementary information, 
an analysis of the correlations between subsequent images allowed confirmation that the variations in a arise 
predominantly from spatial structure rather than dynamical processes within the cell (Supplementary section SI), 
and also allowed the local characteristic length scale of the viscoelastic structure to be determined with 2 nm 
precision. 

Since increased spatial averaging makes small changes in a easier to resolve, an intrinsic compromise is present 
between image contrast and spatial resolution. In this Letter, the choice of 10 nm spatial resolution was found to 
provide sufficient contrast to observe cellular structure. Four example images are shown in Fig. 2, with clear spatial 
structure evident even at the 10 nm spatial resolution limit. The observed spatial dependence shows a number of 
different behaviors because, for each image of a(x), the particle follows a different 3D trajectory. These include 
gradual linear changes in a (e.g. Fig. |2}\) which suggest a spatial gradient in the molecular crowding along the x 
axis, narrow dips in a (e.g. Fig. [2j3 at 40 nm) suggestive of barriers in the cytoplasm, areas of homogeneity (e.g. 
Fig. |2}T), and peaks in a (e.g. Fig.|2j3 at -55 nm) which may follow from small voids in the cytoplasmic structure. 

To characterize the spatial resolution of the imaging procedure and to exclude the possibility that the observed 
spatial variations in a were due to experimental drifts, similar images were constructed for 1 /im radius silica beads 
in water. The resulting images (Fig. |3) show no statistically significant spatial structure, with a = 1 at all spatial 
locations as expected for Brownian motion. Since water is homogeneous with a = 1 throughout, the variation in 
this data allows the statistical uncertainty of our measurements of a to be determined. The spatial resolution was 
characterized by varying the width of the running average along x as shown representatively in Fig. |3]\ to C. As 
the width increases, the spatial resolution is degraded, but the inclusion of more data allows improved precision in 
determining a. The overall precision in a is plotted as a function of resolution in Fig. |3f) for data recorded with 
both squeezed and coherent light, which shows that for fixed sensitivity, the reduced measurement noise achieved 
with squeezing enables 14% enhanced spatial resolution (Fig.[3ji). For instance, from a 10 s data set, biological 
structures that change a by 0.1 are unresolvable using coherent light if their size is smaller than 12 nm. Using 
squeezed light they may be resolved at sizes down to 10 nm. 

Since only the projection of the particle motion onto the x axis is tracked in our current experiments, it is not 
possible to define the complete trajectory along which the ID images are taken. By incorporating our technique 
in a standard 3D particle tracking setup [34| it would be possible to map the complete 3D trajectory. A quantum 
enhanced 3D image of the cell could then be generated at nanoscale, with quantum enhancement only required 



3 



for one axis from which a could be determined. Furthermore, an extended field of view could be achieved by 
increasing the measurement length and thus allowing the particle more time to diffuse, or by applying an optical 
trapping force to the particle to move it into new regions of the cell. 

As is typical of nanoprobe based microscopy techniques, the spatial resolution of the images produced here is 
not fundamentally constrained by the diffraction limit. Since the intracellular medium interacts with the nanopar- 
ticle across its entire surface area, it might be thought that the particle diameter sets a fundamental resolution limit. 
However this is not the case, as can be appreciated from an analogy to atomic force microscopy (AFM). In AFM, 
the transverse resolution is constrained by the tip size. However, the vertical resolution can be far superior[35 1, 
and is fundamentally constrained only by short range forces [ 35 1 and measurement signal-to-noise. Similarly, the 
resolution of the imaging technique introduced here is limited to the particle diameter in directions transverse to its 
trajectory; while in the direction of particle motion far superior resolution is possible, fundamentally limited only 
by measurement noise [27 , 36] and the range of viscoelastic forces from structure in the cell. Consequently, along 
the trajectory of the particle, the location of structure within the cell may be established with precision far smaller 
than the particle diameter. Over time, diffusion acts to drive the particle through the cellular structure on a range 
of trajectories, which combined with 3D tracking, in principle enables nanometer resolution along all three axes. 

The resolution enhancement reported here could be substantially improved using an increased level of squeez- 
ing. Squeezing in excess of 10 dB has been directly measured in a number of experiments Il37ll38l . Using this level 
of measured squeezing, more an order of magnitude enhancement should be feasible, allowing Angstrom level pre- 
cision (Supplementary section S2). In principle, further enhancement may be possible by using an array of photon 
number resolving detectors, with recent theoretical results predicting that this approach could even allow particle 
tracking at the de Broglie limit ||261 . Such a resolution enhancement would be greatly beneficial in biological 
imaging. By comparison, widely used classical super-resolution techniques such as stimulated emission depletion 
(STED) microscopy are currently limited to spatial resolution of a few tens of nanometers within cells ||39ll40l . 

We report the first demonstration of sub-diffraction limited biological quantum imaging, and demonstrate that 
non-classical light can improve spatial resolution in biological applications. The viscoelastic structure within a 
living yeast cell is imaged along the trajectory of a thermally driven nanoparticle, revealing spatial structure with 
length scales down to 10 nm. The use of squeezed light was found to enhance the spatial resolution by 14%. 
Future experiments which apply this method with improved squeezing technology are expected to resolve sub-nm 
structure. 
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Figure 1 : Experimental setup. (A) Counter-propagating trapping fields (orange) confine particles between two ob- 
jectives, and are isolated from the detector with polarizing beamsplitters (PBS) and waveplates (A/2). An imaging 
field (green) allows visual identification of the particles near the optical trap on a CCD camera. The particle track- 
ing measurement relies only on an amplitude squeezed local oscillator and an amplitude modulated probe (red), 
with the probe providing dark-field illumination, and the particle tracking signal arising from interference between 
scattered light from the probe and the local oscillator. (B) Measured particle motion, which is the x projection of 
the 3D motion (shown schematically in the inset). (C) The MSD is constructed with both squeezed light (dark red) 
and coherent light (gold), and a determined by fitting this to Eq. [2] The classical and squeezed example traces 
here both yield a = 0.83. The normalized spectral density is plotted in the inset, which shows that squeezing 
suppressed the noise floor by 2.4 dB. (D) The raw data was divided into 100 ms segments and the value of a 
established for each (solid dots). The light red shaded region represents the moving mean and standard error with 
a 0.5 second width. 
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Figure 2: Images of a. Each circle represents a single measurement of a vs x using a 100 ms set of data. The 
shaded regions represent the running mean and standard error with 10 nm resolution (thick black bar). The images 
(A-D) were recorded minutes apart, which allowed the particle time to diffuse to different regions of the cell, with 
qualitatively different spatial structures. The particle confinement is greatest where a is lowest, such as the dip 
about 40 nm in (B), and the particle movement is most free when a is highest, such as the peak at -55 nm in (D). 
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Figure 3: To calibrate the quantum enhancement achieved here, an image of a was constructed from particles 
in water. The measurements of a vs. x are shown in (A). (B-D) show the corresponding moving mean and 
standard error calculated over a 200 nm range with spatial resolution of 2 nm, 10 nm, and 50 nm, respectively. The 
data closely follows the expected a = 1 result (horizontal line). The number of independent spatially resolved 
points decreases as the spatial resolution is broadened, while the increased spatial averaging allows a more precise 
determination of a. (E) The precision with which a could be determined is plotted as a function of spatial resolution 
for data recorded with both squeezed (solid line) and coherent (dashed line) light. Since the number of points 
being averaged is proportional to the spatial resolution, the precision scales as the inverse square root of spatial 
resolution until the averaging window width becomes comparable to the spatial range of the measured data. When 
resolving features of a fixed amplitude in a, squeezing allows enhanced spatial resolution as shown in (F), with 
approximately 14% narrower spatial resolution available over most of the measurable range. Outside the vertical 
dashed lines in (A-F), there was insufficient data to construct an image. 
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Supplementary Materials 

SI: Correlation analysis of images 

It is important to rigorously establish whether the changes in a measured here result from particle motion within 
nanoscale structure or from time-dependent cellular processes. Because the trajectory over which a is imaged 
differs in each run of the experiments reported here, different images correspond to different regions of the cell. 
However, when spatial images are produced within a short time of one another, the corresponding trajectories are 
not spatially separated by a substantial distance, and the images should retain some correlations. By contrast, 
any temporal contributions to the measured a(x) should be completely uncorrected between subsequent images. 
Thus, analysis of the correlations between subsequent images can determine whether temporal fluctuations or 
spatial structure in a are dominant. The correlations which result from spatial structure can be predicted with 
a simple model which assumes that a is perfectly correlated if the displacement dr along the y and z axes is 
below the length scale of structural variations r c , and completely uncorrected if dr > r c . Then, assuming that 
the motion along each axis independently follows Eq. [2 and also assuming Gaussian probability distributions, the 
average correlation is given by the probability that the particle remains within r c 



where T c is the average time that the particle takes to diffuse into an uncorrected region of the cell, such that 
(dr 2 ) = rj?. To determine whether the experimental data follows Eq. [3] a series of a(x) images were constructed, 
and the correlations between the subsequent images calculated as a function of delay (Fig. SI). This experimental 
data was fitted to Eq. [3] with T c and a used as fitting parameters. The data shows excellent agreement with the 
simple model, which confirms that the measured changes in a are spatial rather than temporal. This is consis- 
tent with observations from other types of cell that have found the viscoelasticity to vary spatially by an order 
of magnitude lETTl . which is much larger than the temporal variations produced through cellular processes lUTl . 
Furthermore, by quantifying the time constant T c , the characteristic local length scale of viscoelastic structure can 
be determined in the region of the nanoparticle. To test whether this length scale varies across the cell, the corre- 
lation analysis described above was carried out with three separate sets of data, as shown in Fig. SI. Assuming 
the diffusion follows the MSD measured in Fig. [2f , the length scales of the viscoelastic structure in each region 
were determined to be 46.9 ± 1.3 nm, 43.7 ± 2.4 nm, and 42.6 ± 1.9 nm. This demonstrates that changes in the 
characteristic length of spatial structure in different parts of the cell can be statistically distinguished. 

S2: Theoretically predicted quantum enhancement 

Here the enhancement in spatial resolution achievable with squeezed light is determined by analyzing simulated 
data. The data analysis in the main text is performed as follows. A discrete series of measurements x n of the 
particle position x are performed at evenly spaced intervals dt in time, where n is the measurement number. These 
are used to determine the mean squared displacement (MSD), by the relationship 



where -/Vmsd ls the mean variance. Although the expectation value of this experimental MSD is given by 2Dr a 
(Eq. 2 of the main text), a single trajectory will only asymptote towards this in the limit of infinite statistically in- 
dependent measurements. The number of statistically independent measurements at any delay r is given by T/t, 
where T is the total measurement time. As such, even a perfectly noise-free determination of the particle motion 
cannot yield a with complete certainty in a finite measurement time. Thus, the stochastic nature of the MSD intro- 
duces an intrinsic uncertainty into the value of a in addition to any uncertainty associated with the measurement 
noise. This intrinsic uncertainty is greatest at long delays, where the number of independent measurements is 
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Figure SI: Correlation between a series of measured sets of a(x) as a function of the time separation between the 
data sets. The correlation decays in time because of the unknown motion along the y and z axes. The circles are 
experimentally determined correlations between a series of measurements, and the lines fit the data to Eq.|3] which 
determines the probability that subsequent measurements sample the same spatial environment, a was optimized 
to minimize the mean square error for a range of T c . From this the optimal T c was found, with the uncertainty 
corresponding to the change in T c which doubles the mean square error. 
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Figure S2: The influence of measurement noise in viscoelastic imaging. (A) The precision with which a can be 
determined {5a) is plotted against the measurement noise variance, normalized to the MSD at minimum delay. 
The dotted line scales as the square root of the noise variance, and approximately follows 8a between variances 
of 10 and 100. The inset zooms in on this result in the regime of the experiments here, with the vertical dashed 
lines showing the noise variance on the experimental data. (B) The ratio of the spatial resolutions achievable with 
squeezed light and coherent light, with either 2.4 dB (blue) or 10 dB (green) of squeezing. 
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lowest. Measurement noise, by contrast, is of approximately constant amplitude with delay, and therefore has the 
greatest relative effect at short delays where the expectation value of the MSD is smallest. 

The precision of any determination of a depends on both the intrinsic stochastic uncertainty and the measure- 
ment noise, such that the enhancement in a precision which may be achieved via squeezed light cannot be simply 
related to the noise suppression. To quantitatively predict the enhancement achievable, we perform Monte Carlo 
simulations, in which 10000 particle trajectories are simulated by a normally distributed random walk. Varying 
levels of white measurement noise are added to the trajectories, after which the MSD is calculated and a deter- 
mined. The uncertainty da is then given by the deviations of the predicted a from the ideal result of a = 1 which 
should follow from such a random walk. 

By repeating this for various levels of measurement noise, we could plot the uncertainty in a as a function of 
the measurement noise variance (Fig. S2A). In this case, the measurement length and fitting range matched the 
measurement conditions of the main text, and the noise variance was normalized to the MSD at minimum delay. 
When the noise is far smaller than the MSD (below 1 in Fig. S2A), the stochastic walk dominates the uncertainty 
in a, and techniques for noise reduction such as squeezing yield little improvement. When the measurement noise 
dominates the MSD at all analyzed delays, a cannot be determined with any reliability. As such, this region is not 
included in in Fig. S2. Between these two extremes, measurement noise only dominates the MSD at short delays, 
and the precision with which a can be determined scales approximately as the square root of the noise variance. 
The relative contribution of shot noise increases as the measurement rate increases or the measurement duration 
decreases, both of which are required to produce high precision images, or if the optical power is reduced. As 
such, quantum enhanced precision provides maximal enhancement in the regime which is both most challenging 
to reach in classical experiments and most relevant to biological imaging. 

When performing spatially resolved experiments, the precision with which a can be determined at a specific 
location is given by Sai magc — SaN^ 1 / 2 , where N is the mean number of data points averaged at each location. 
For ID imaging, this is given by 

N = N tot ^, (5) 
a^FOV 

where N to t is the total number of measurements of a, xfov is tne field of view, and dx is the spatial resolution. 
Thus, if a fixed precision <5ai mago is required in a ID image, the spatial resolution scales as dx oc 5a 2 . Using this, 
the spatial resolution enhancement is plotted in Fig. S2B for both 2.4 dB and 10 dB of measured squeezing. 

At the noise levels measured here, we can predict that 2.4 dB of measured squeezing should have allowed 
21% narrower spatial resolution in our experiment, which is greater than the measured 14% enhancement. The 
discrepancy between predicted and achieved enhancements may be due to the technical noise in the experiment 
such as low frequency noise and unwanted inclusion of RF modulations that were used for phase locking of the 
squeezing cavity. If the measurement rate were increased, such that the shot noise variance was 100 times greater 
than the MSD at minimum delay, 2.4 dB of measured squeezing would allow a 45% reduction in spatial resolution. 
If state-of-the art squeezing were used and the measured squeezing reached 10 dB, the quantum enhancement could 
then provide 90% narrower spatial resolution. Alternatively, an even greater measurement rate could sufficiently 
increase the measurement noise to preclude determination of a from a classical experiment, while a could still be 
determined with squeezed light. 
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